Charge dynamics and Kondo effect in single electron traps in field effect transistors 
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We study magneto-electric properties of single electron traps in metal-oxide-semiconductor field 
effect transistors. Using a microscopic description of the system based on the single-site Anderson- 
Holstein model, we derive an effective low energy action for the system. The behavior of the system 
is characterized by simultaneous polaron tunneling (corresponding to the charging and discharging 
of the trap) and Kondo screening of the trap spin in the singly occupied state. Hence, the obtained 
state of the system is a hybrid between the Kondo regime, typically associated with single electron 
occupancy, and the mixed valence regime, associated with large charge fluctuations. In the presence 
of a strong magnetic fleld, we demonstrate that the system is equivalent to a two level-level system 
coupled to an Ohmic bath, with a bias controlled by the applied magnetic field. Due to the Kondo 
screening, the effect of the magnetic field is significantly suppressed in the singly occupied state. We 
claim that this suppression can be responsible for the experimentally observed anomalous magnetic 
field dependence of the average trap occupancy in Si — Si02 field effect transistors. 

PACS Numbers: 73.20.-r, 71.38.-k, 75.20.Hr 



I. INTRODUCTION 

Experimental techniques probing dynamics of a few 
state quantum system have been of great interest re- 
cently as they may provide new prospects for the devel- 
opment of electronics and computing. Examples of such 
systems and techniques in solid state include quantum 
dots^, superconducting qubits^, magnetic resonant force 
microscopy'^, and single electron traps in a metal-oxide- 
semiconductor field effect transistor (MOSFET)"*"*. 

This latter system, which is a subject of this work, con- 
sists of a defect located near the oxide-semiconductor in- 
terface of a MOSFET (Fig 1(a)). The tunnehng between 
the defect and the 2-dimensional electron gas (2DEG) 
in MOSFET inversion layer manifests itself as a random 
telegraph signal (RTS) in the transport current through 
the MOSFET. When a certain trap energy level (whose 
energy can be controlled by the gate voltage) crosses the 
chemical potential of the 2DEG, electrons can hop be- 
tween the level and the conduction channel, thus charg- 
ing and discharging the trap. For a sufhciently small 
MOSFET's this leads to the sudden switching in the re- 
sistivity of the conduction channel and hence to RTS in 
the transport current. 

The charge dynamics of the trap exhibits a number of 
features consistent with dynamics of a two-level system 
(TLS) coupled to an Ohmic environment^. In particular, 
the experimental dependence of tunnel rates on the TLS 
bias (controlled by the gate voltage) and external tem- 
perature was found to agree with those calculated for the 
spin-boson Hamiltonian^. However there is a number of 
RTS properties that cannot be explained based on the 
TLS phenomenology. 

First of all, the RTS is observed to occur on a relatively 
large, millisecond-to- second tivcie scale, which seems to in- 
dicate that the observed traps are positioned sufficiently 
far, at distances 20-25 A from the 2DEG. At the same 
time, the direct electrostatic measurements of trap posi- 



tions indicate that this is usually not the case - the traps 
are located only a few A from the inversion layer*. 

Another aspect that clearly lies outside the scope of 
TLS phenomenology has recently been revealed in the 
experiments by Xiao et af , where the dependence of 
the RTS on applied magnetic field was studied. In 
these experiments, from the gate voltage and the mag- 
netic field dependence of RTS, it was determined that 
the most likely origin of RTS was the random switch- 
ing between (spin-1/2, neutral) and (spin-0, negatively- 
charged) states of the trap. In particular, at sufficiently 
high temperatures (above several degrees K), the proba- 
bility of the filled state of the trap was rapidly decreas- 
ing with applied magnetic field, consistent with the sim- 
ple model of a spin-1/2 empty state and a singlet filled 
state. However, at lower temperatures significant devi- 
ations from the simple paramagnetic behavior was ob- 
served, possibly indicating quenching of trap's magnetic 
moment''. Such magnetic behavior appears to be consis- 
tent with the Kondo effect. 
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FIG. 1. (a) Schematics of traps in MOSFETs; (b) A dia- 
gram for the trap energy levels. Coupling to optical phonons 
shifts the twice occupied level. 
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However, it is well known that the Kondo screening is 
only effective for siifficiently strong hybridization of the 
localized state with the continuum. Given the extremely 
slow observed tunnel rate, one should conclude that the 
Kondo temperature should be negligibly small, thus rul- 
ing out the possibility of the Kondo effect explanation. 

The purpose of this work is to show that the slow 
charge dynamics and the Kondo effect are in fact not 
mutually exclusive, if we take into account the strong 
electron-optical phonon coupling in the oxide layer of a 
Si MOSFET. In our previous work^*^ we have demon- 
strated that this coupling is the origin of exponential 
renormalization in RTS timescales and explains the in- 
consistency between expected and observed positions of 
traps. In the present work, starting from the microscopic 
description of the system based on Anderson-Holstein 
model, we derive an effective low energy action for the 
charge and spin dynamics of the trap. We find that the 
low-temperature behavior of the system is characterized 
by simultaneous polaron tunneling (corresponding to the 
charging and discharging of the trap) and Kondo screen- 
ing of the trap spin in the singly occupied state. In the 
presence of a strong magnetic field, we map our system 
onto a two level-level system (TLS) strongly coupled to 
an Ohmic bath, with a bias controlled by the applied 
magnetic field. Due to the Kondo screening, the bias in- 
troduced by the magnetic field is significantly suppressed 
compared to the bare Zeeman energy. We claim that this 
suppression can be responsible for the experimentally ob- 
served'' anomalous magnetic field dependence of the av- 
erage trap occupancy in Si — Si02 field effect transistors. 
It is crucial that the Kondo temperature need not be 
small as it is controlled by the bare hybridization matrix 
element between the localized and conduction electrons, 
i.e., it is not affected by the polaronic slow-down. 



II. MODEL 

To describe the electronic part of the trap-channel sys- 
tem we use the Anderson Hamiltonian: 

Ha = ^Ekcl^Cka- + ^Eodldo- + Un^ni 

(7 

+ Y,T,a{clj, + dlck,), (1) 



k a 



k a 

where cl,^{cka) and dl,{dcr) are creation(annihilation) op- 
erators of electrons in the conduction band and at the lo- 
calized defect orbitals respectively, rio- = rf^d^ =t, i)- 
U and Ted are the on-site Coulomb energy and the hy- 
bridization matrix element respectively. The coupling 
strength between the localized states and the 2DEG can 
be characterized by F = TrvT^^, where v is the density 
of states in the 2DEG. The single particle level of the 
trap is assumed to be positioned deep below the chemi- 
cal potential of 2DEG, Eq < fi, as shown in Fig. 1(b). 



For simplicity we will assume that the conduction band 
is symmetric and set /i = 0. 

Another interaction that we include in our model is due 
to optical phonons in the Si MOSFET oxide layer. Si02 
is a polar crystal and known to exhibit strong electron- 
optical phonon coupling leading to formation of polaronic 
states. We incorporate this electron-phonon coupling in 
our model as 



2 2^2 
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Here x and p are displacement and conjugate momentum 
of local optical phonon of frequency loq {loq ~ 50 meV at 
Si — Si02 interface) , m is the phonon effective mass (of 
the order of Si02 crystal unit cell mass), and A is the 
coupling constant between the excess charge in the trap 
and the phonon. We assume the trap state with one 
electron is neutral and therefore introduce the off-setting 
— 1 term in the electron-phonon interaction in Eq. (2). 
The strength of the electron-phonon coupling has been 
estimated in Ref.^° and can be expressed in terms of the 
polaronic shift Ep = j^muj^ ^ 1 eV. In the following 
we will assume that U > Ep ^ T > luq. 



III. EFFECTIVE ACTION 

We analyze low temperature partition function of the 
system Z = Tr[exp— where H is given by Eq. (2) 
and l3~^ is smaller than any energy scale in the sys- 
tem. It is convenient to decouple the U term in Hamil- 
tonian (1) by means of Hubbard-Stratanovich transfor- 
mation by writing the spin-spin interaction as Uni^ni = 
-{U/2){n-(-nif + {U/2){n^+niy^. The partition func- 
tion can be then cast in the functional integral form as 
Z = 



I 
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[X{r)-Y(r)] 



), (3) 



where 7?.[^(r)] = Z(r) 4(T)d„(T), A = Eq + (C//2) + 
2Ep. Here the explicit "time" dependence of dJ.(do-) is 
defined by the interaction representation with respect to 
kinetic + tunneling terms in the Hamiltonian (1). In 
Eq. (3) we have introduced X = Xx+Eo+{U/2) as a coor- 
dinate of the oscillator in units of energy and M = m/X^ 
as the oscillator's mass in the corresponding imits. The 
scalar field Y is conjugate to the spin of the trap, T 
stands for time-ordering and the brackets (•) denote aver- 
aging with respect to electronic degrees of freedom. The 
averaging of the ordered exponents in Eq. (3) can be done 
by following the prescription of Refs.^^'^'^ based on the 
theory of singular integral equations. Here we provide 
the result: the functional integral in Eq. (3) can be cast 
in the form Z = J V[X, Y] exp [- dT{V + K.)], where 
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functionals V and K. represent potential and kinetic en- 
ergy of a particle moving in two-dimensional (X,Y) plane: 



V{Z) = I - ^tan- 



27r 



log 



(4a) 



K = ^— + K{X + Y) + K{X - Y) , 



Z{t) 

2-2(t) - Z2(t') °^r2+Z2(T 

The potential energy of the "particle" is presented in 
Fig. 2(a). It has four local minima. The two minima on 
the X axis (at X2 = A-AEp and Xq = A for U, Ep » T) 
correspond to different charge occupations of the trap, 
i.e., by 2 and by electrons. 



dr 

r2 + Z2(t) 
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FIG. 2. (a) Local minima for the effective potential in 
Eq. (4b); (b) A classical path contributing to the transition 
amplitude between a 2 electron and a spin-up electron states; 
(c) Example of a path containing a spin-kink "monopolc". 
Such trajectories do not contribute to the transition ampli- 
tudes; (d) Trajectory that involves multiple charge-kinks 



Up to an additive constant (which we have dropped 
in Eq. (3)) potential energies of these minima are V2 = 
2Eo + U — Ep and Vo — ~Ep. The two minima off the X 
axis (at X^ = A - 2Ep, = (7 and = A - 2Ep, 
Yl = —U) are degenerate in the absence of external 
magnetic field with Vf = Vj = Eq and obviously cor- 
respond to the occupation of the trap by 1 electron with 
up and down spin, respectively. As can be seen directly 
from Hamiltonians (f ,2), the minima can be interpreted 
as zero temperature free energies of the trap occupied 
by 2, and by 1 electrons in the limit of vanishingly 
small F. In the RTS experiments the traps are filled 
with either 1 or 2 electrons^, and therefore we assume 
that Vt =i ~ V2 < Vo (or Eoc^2Eo + U - Ep < -Ep). 



IV. EVALUATION OF FUNCTIONAL 

The kinetic part of the energy in Eq. (4) will force the 
"particle" to tunnel between the minima. Because the 
energy of the empty state (denoted by the blank circle in 
Fig. 2(a)) is greater than that for the other three states 
(the shaded minima in Fig. 2), it will be effectively de- 
coupled from those states. (Transition amplitudes to the 
empty state will contain exponentially small prefactors 
exp — /3(Vo — V2) <Si f.) An example of the "classical" 
trajectory for a transition between the two electron state 
and spin-up state is presented in Fig. 2(b). Let us first 
write down explicitly the transition amplitude that cor- 
responds to the first diagram in Fig. 2(b). In doing so 
we must, in principle, determine the classical trajectory 
for the effective action given by Eq. (4) between the two 
local minima corresponding to the singlet and spin-up 
states. Solving the resulting classical equation of mo- 
tion, however, is difficult. Instead we approximate the 
classical path by a piece-wise linear kink trajectory, as 
shown in Fig. 2(b)^2,i3^ Duration of each "charge" kink 
(i.e., transition between the 2 electron state and a 1 elec- 
tron state), which we denote by t'^, can be determined 
by minimizing the effective action for a single linear kink 
trajectory with respect to t'^. Upon substitution of such 
linearized trajectory in Eq. (4) the transition amplitude 
of interest can be written as 



Jo JO ^ 



ji5c(-rf-r|)-a;c log ■ 



where parameters of and r'^ are 

(2 _.U^^ 
ac = - tan — ■ 
\7r i 



(5) 

(6a) 
(6b) 
(6c) 



In the above expression for the transition amplitude the 
kinks at time moments rf and T2 interact logarithmically 
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and Qfc is dimensionless coupling strength. The fugacity 
which is essentially the probability for a kink to oc- 
cur at a given time, is suppressed by an exponentially 
small quantity in Eq. (6b) as a consequence of the os- 
cillator having a finite mass. This suppression of charge 
tunneling is in agreement with the our previous result 
where we argued that the inconsistency between the ob- 
served and expected tunnel rates in the RTS systems 
is due to the presence of a local phonon strongly cou- 
pled to the trap charge. The quantity 5c in the above 
equation is a bias between one- and two-electron states, 
<5c = V2 - Vta =Eo + U -Ep. 

Next let us look at the transition amplitude that in- 
volves spin transitions. In Fig. 2(b) this amplitude can 
be singled out as an effective self-energy part of the dia- 
grams that involve 2n, (n = 0, 1, ...) kinks and anti-kinks. 
Applying the same procedure, i.e., assuming that all the 
trajectories are piece-wise linear, evaluating their contri- 
bution to the effective action in Eq. (4) and minimizing 
the action with respect to the width of the spin kinks , 
we obtain that the spin transition amplitude, -Z^ in Fig. 
2(b), is given by 

4 / — ^- / -T / 



X exp 



X exp 
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(7) 

(8a) 

(8b) 
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Here, bg = QI^-bE is the energy splitting between the bare 
spin-down and spin-up states due to the applied mag- 
netic field B^^. The transition amplitude Z'^ in Eq. (7) 
is of the same form as the grand partition function of 
a Coulomb gas of ± kinks positioned along a straight 
line of length {r^ — rf ) interacting logarithmically with 
the coupling strength and chemical potential defined 
as e"'* = In order to evaluate the full transition 

amplitude in Fig. 2(b) we must, in principle, take into 
account that the charge kinks at the ends of the inter- 
val, i. e., at "positions" and r^, interact (also loga- 
rithmically, as can be seen from the effective action in 
Eq. (4)) with the spin kinks at positions rj*. One can 
see, however, that this interaction is effectively weak. In- 
deed, kinks at rj* are coupled strongly with each other 
and tend to form close pairs (dipoles) of average effec- 

tV(2 - OLs) ~ 1/r, while 



the dipoles are separated, on average, by distances of 
the order of I ~ rV^^ _ [//p2 15^ r^^^^ YwlVh at 
are separated from those at r| by a distance (time) at 
least of the order t'^ ~ l/t^o, which physically repre- 
sents the time needed for the formation of the dressed 
electron-phonon state in the trap. Thus, the interac- 
tion energy for a kink at rf with a dipole at rf -|- r is 
of the order rf/r. Summing over the dipoles (with the 
closest dipole located at distance of the order from 
the charge kink at and the farthest dipole roughly 
at rl) one finds that the charge kink-dipole interac- 
tion energy for the transition amplitude in Fig. 2(b) is 

~ {djl) j;r^^ dr/r ~ (r/[/) log [(r| - rf )/Te]. This in- 
teraction is of the same long-range form as that between 
the kinks at and t^. e.g.. Eq. (5), however paramet- 
rically it is much weaker, ofg ~ 1 ^ (T/U). Therefore, 
interaction between the charge kinks and the spin kinks 
can be neglected. Also, we do not need to consider con- 
figurations with an odd number of spin kinks between the 
charge kinks, e.g. Fig 2(c). That is because the energy of 
an unpaired kink diverges logarithmically with (r^ — rf ) , 
and hence its formation is not favorable. Thus the tran- 
sition amplitude that contains two charge kinks at rf 2 
can be written as 



where Zk 



dr^ 



-■K- 



We are now ready to extend this expression to a higher 
order amplitudes. One such amplitude with four charge 
kinks at is shown in Fig. 2(d). We notice that only 
charge kinks connecting the singlet state with the same 
spin state interact. Hence, wc; assign superscript a to 
distinguish between different types of charge kinks, rj —^ 
Tp . Using the same energetic arguments as before, we 
can show that the interaction between the spin kinks in 
different domains, as well as the interaction between any 
charge kinks and spin kinks, is negligible. Therefore, the 
transition amplitude can be written as 

n 



X exp 
X exp 



2n 



-Sc^i-iyr; 

i=i 

2n _ 



(9) 



In zero magnetic field (Sg = 0), the spin partition func- 
tion zj^(T - r') = zj^(r-r') ~exp [-(r - r')-Fo], where 
Fq is the free energy per unit length for the Coulomb gas 
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described by the classical partition function of Eq. (7). 
In this regime, the full partition function of Eq. (9) corre- 
sponds to a gas of charges of two flavors, with charges at 
points T2j and T2j-i having the same flavor and opposite 
charge, and only charges of the same flavor interacting 
with each other. In the presence of applied magnetic 
field, the symmetry between the two fiavors is broken. 



V. MAGNETIC FIELD DEPENDENCE 

To evaluate the effect of magnetic field on partition 
function Eq. (9), we first determine the dependence of 
Z]^ on magnetic field. For that we use the scaling proce- 
dure for the Coulomb gas due to Anderson et af^. The 
renormalization group equations read: 



rflog 



-ia,e.+0{(t) , 



d log T* 



(10a) 
(10b) 
(10c) 



ZK[i/3/Tn,^s,a,]=e^p{pAF)ZK[ifi/fn,^s,a,] . (lOd) 

where AF = J^^ dlogr^'^g/r*'. For sufficiently strong 
magnetic field, {5s > Tk = t*~^ exp [— 1/^^]), the scaling 
should bo terminated when 5gf^ ^ 1. At this point the 
magnetic field becomes very strong (on the renormalized 
energy scale), while the long-range coupling constant 
decreases and so the system moves away from the quan- 
tum phase transition point at a., = 2 into the disordered 
phase^^. For mathematical simplicity we will assume that 
2 — as— 4^s in Eqs. (8). In such symmetric case the scal- 
ing Eqs. (10b) and (10c) coincide (up to quadratic terms 
in and 2 — as) and can be easily solved analytically. 
For <C 1 the resulting renormalized Sg and are 

~5s = 5s cxp[-(l/2) \og-\5s/TK) + 0(6)] , (Ha) 
I = log-i {5s/Tk) + 0[log-3 {5s/Tk)] . (lib) 

It is easy to see that in the rescaled system the long 
range coupling becomes irrelevant. A simple perturbative 
estimate for the average size of a dipole is f */(5s'f*)^ ~ 
and therefore both intra- and inter-dipole interactions are 
negligible. Zk can then be readily summed by means of 
Laplace transform^° with a simple result 

Z^(r)= exp{r[F^ + ~5s + 0{\og-^ {Ss/Tk))]} , 
^i,(T)« es exp {t[F^ + 5s]} + exp {t[F^ -ds]}, (12) 

where Fq is magnetic-field-independent free energy and 
5s is given by Eq. (11). In the long time limit that 
we are interested in, is negligible compared to Zj^. 
Hence, in the partition function Eq. 9, only the kinks 
that correspond to the transition between the singlet 



and the spin-up states survive. The partition function 
becomes then equivalent to the one for a two-level sys- 
tem coupled to an Ohmic bath®. Eq. (12) together with 
Eq. (9) indicate that an external magnetic field intro- 
duces an effective bias for the TLS. This bias, how- 
ever, is significantly reduced (see Eq. (11)) as compared 
to the bare Zecman energy of a paramagnetic spin due 
to the strong coupling between kinks in Eq. (7). The 
partition function of the latter model is equivalent to 
that of the Kondo model^^ and, therefore, the suppres- 
sion of the magnetically-induced bias of the TLS is es- 
sentially the Kondo effect. Indeed, the magnetization 
in the singly occupied state can be easily computed as 
M « d5s/d{giiBB) = i(l - l/[21og(ff/xsB/TK)]), con- 
sistent with the Bethe ansatz result of Andrei et al}^. 



VI. DISCUSSION AND CONCLUSIONS 

It is important that the Kondo energy scale Tk is not 
determined by the extremely slow charge tunneling rate 
(of the order of £,c/t'^)- Rather, the fugacity ^s, which 
determines Tk in Eq. (12), is due to virtual transitions 
between Fermi sea and the localized orbital. In the lan- 
guage of the effective action in Eq. (4) these virtual tran- 
sitions correspond to timneling of the y-field, which, un- 
like the X-field, is massless. Therefore, as can be seen 
from Eq. (7), the fugacity does not contain the ex- 
ponentially small quasiclassical suppression factor corre- 
sponding to the "penetration under the barrier" , which 
strongly reduces the charge fugacity (Eq.(6)). As a 
result the Kondo temperature (~ exp[— l/^^]) can be 
quite large leading to significant renormalization of the 
magnetic-ficld-induced TLS bias. 

Based on the above presented arguments, we believe 
that the anomalously weak magnetic field dependence of 
the interface trap occupancy observed in the low tem- 
perature experiments of Xiao et al. can be explained by 
the Kondo screening of the local moment in the singly 
occupied state of the trap, assuming that the Kondo 
temperature is Tk ~ 1 K. Then, for low temperatures, 
T < Tk, we expect that the ratio of the probabilities cor- 
responding to single (spin-up) and double occupancies of 
the trap should scale as P1/P2 oc exp{giJ.BB/2T), where 
the effective magnetic field B is related to the applied 
magnetic field B as B = Bexp{-l/[2log{giJ.BB/TK)]) 
for gusB > Tk- 

Clearly, further detailed experimental study is neces- 
sary to test out theoretical picture for the Kondo effect 
in the electrically active defects (traps) in Si field ef- 
fect transistors. It is worth noting, however, that similar 
physics, that is simultaneous slow charge dynamics and 
the Kondo effect, is expected to occur in any other sys- 
tem that has defects located at the interface between a 
strongly polar insulator and a conductor. One possible 
way to detect the Kondo effect is to look for low temper- 
ature anomalies in the FET channel resistivity. Another 
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approach is to look for spectral features (Kondo peak) by 
direct tunneling through a defect at the Si02-Si interface. 

In summary, we have derived an effective action for the 
charge dynamics in a single electron trap in a Si MOS- 
FET and have shown that in the limit of a strong mag- 
netic field [gfiBB > Tk) it is equivalent to a two level 
system strongly coupled to an Ohmic environment. The 
effective bias of the TLS can be controlled by an external 
magnetic field. However, the magnetic field dependence 
of the trap spin at low temperatures is suppressed due to 
the Kondo effect. 
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